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Abstract 

Recently various new concepts for the construction of Dirac operators in lattice Quantum Chromodynamics (QCD) 
have been introduced. These operators satisfy the so-called Ginsparg- Wilson condition (GWC), thus obeying the 
Atiyah-Singer index theorem and violating chiral symmetry only in a modest and local form. Here we present 
studies in 4-d for SU(3) gauge configurations with non-trivial topological content. We study the flow of eigenvalues 
and we compare the numerical stability and efficiency of a recently suggested chirally improved operator with that 
of others in this respect. 
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1. Motivation and technicalities 

A priori continuum Dirac operators have chiral 
symmetry, which then is broken spontaneously in 
the full QCD dynamics. Lattice discretizations of 
the Dirac operator are bound to violate the chiral 
symmetry. However, if they obey the GWC this vi- 
olation is in its mildest form. Currently three types 
of exact realizations of so-called Ginsparg- Wilson 
(GW) operators are known: The overlap opera- 
tor, perfect actions and domain wall fcrmions; for 
reviews on recent developments cf. Q. In Q we 
suggested a method to systematically expand the 
lattice Dirac operator in terms of a series of sim- 
ple lattice operators and to turn the GWC into 
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a large algebraic system of coupled equations for 
the expansion coefficients. The solution for a fi- 
nite parametrization leads to an approximate GW- 
operator. Such a chirally improved Dirac opera- 
tor has been constructed in 2-d for the Schwinger 
model and also for 4-d QCD , where an ana- 
lysis of its spectrum revealed intriguing properties. 

The spectral properties of QCD Dirac opera- 
tors provide an efficient means to study proper- 
ties of the QCD vacuum related to chiral sym- 
metry breaking. The spectral density in the limit 
of vanishing eigenvalues is proportional to the or- 
der parameter of chiral symmetry breaking via the 
Banks-Casher relation Q . The microscopic distri- 
bution density exhibits features in a universality 
class discussed in random matrix theory The 
eigenmodes for exactly vanishing eigenvalues ap- 
pear to represent instantons (classical solutions of 
the QCD field equations). In the quantization of 
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a lattice-regularized field theory it may come as a 
surprise that instanton-like objects (originally de- 
fined on differentiable manifolds) may play a role. 
Indeed, they are a miniscule contribution to the 
total action of a gauge configuration. However, as 
will be seen, they are observed in realistic gauge 
field configurations and they characterize the sec- 
tor of small eigenmodes R |j . 

The overlap operator H has the form 



Dov = I — Z with Z = 75 sign(i7) , 



(1) 



where H is related to the hermitian Dirac oper- 
ator. This is constructed from an arbitrary Dirac 
operator Dq (like e.g. the Wilson operator), H = 
75 (s — -Do); and s is a parameter which may be 
adjusted in order to minimize the probability for 
zero modes of H. In case Dq is already an overlap 
operator one reproduces Dq^ — Dq for s = I due 
to sign (sign(i/)) = sign(i/). 

The sign-function may be defined through the 
spectral representation, but this is technically im- 
possible for realistic QCD Dirac operators, which 
for lattices have dimension 0(10''^ — 10^) (~ 
ncoior'noirac' L^)- Notc, that in the subsequent ap- 
plications the operator has to be determined many 
times, since it may be entering a diagonalization 
or a conjugate gradient inversion tool. One there- 
fore relies on the relation sign{H) ^ H/ \J and 
approximates the inverse square by some method. 
In our computations ||lO| we follow the methods 
discussed by |ll|, 12, j^. One approximates the 
inverse square root by a Chebychev polynomial, 
which has exponential convergence in [e, 1], where 
e (and thus the order of the polynomial) depends 
on the ratio of smallest to largest eigenvalue of 
H^. Clenshaw's recurrence formula provides fur- 
ther numerical stability. 

Depending on the input Dirac operator Dq the 
rate of convergence may be quite unfavorable due 
to small eigenvalues of . Therefore one computes 
the inverse square separately for the subspace of 
the smallest 20 eigenvalues (using the spectral rep- 
resentation) and the reduced operator (with poly- 
nomial approximation) . The subspace has to be de- 
termined with high accuracy. All diagonalizations 
have been done with the Arnoldi method H . 




Fig. 1. We display the flow of eigenvalues of a matrix Z 
derived from a randomly constructed matrix H with all 
eigenvalues ±1 except for one, which is slowly changed 
between and —1 (begin and end positions: full circles, 
value 0: square, intermediate values: dots). 



2. Flow of eigenvalues 

Lattice Dirac operators are necessarily non- 
hermitian, but "75-hermitian" , i.e. one may always 
define a "hermitian Dirac operator" 75 D. The 
eigenvalues of a realistic "hermitian (e.g. Wilson-) 
Dirac operator" are distributed broadly on the 
real axis and have also values close to the origin. 
The sign-operator in (1) projects all n+ positive 
eigenvalues to -1-1, and all n_ negative eigenvalues 
to —1. The difference — n_ is twice the number 
of zero modes of Dov 

We are interested in the case when (1 — Dq) is 
close to, but not quite an overlap operator. In par- 
ticular we want to study here the emergence or 
vanishing of zero modes. As a test case we con- 
struct a random, hermitian matrix operator like 
H with eigenvalues ±1. For 75 we choose an ar- 
bitrary realization like = (~1)* ^ij- The eigen- 
values of the corresponding Z are then distributed 
randomly on a unit circle around the origin. Due 
to "75-hermiticity" of Z its eigenvalues occur in 
complex conjugate pairs. Real eigenvalues are at 
positions ±1 giving rise to vanishing eigenvalues of 
Dov, hence we call them zero modes. 

Starting with a matrix where ri+ = n_ we will 
find no such "zero modes" . Let us change contin- 
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Fig. 2. Real eigenvalue x associated with the instanton 
eigenmode vs. instanton radius in lattice units: Wilson 
operator, overlap operator (with Dq the Wilson operator), 
and the chirally improved operator. 

uously one of the eigenvalues from its value A = 
+1 to —1 (Fig.l). Starting with all eigenvalues on 
the unit circle (none of them at the real positions 
±1) a complex conjugate pair moves to the inte- 
rior. When A changes sign, the pair becomes a pair 
of two real eigenvalues, moving towards the circle 
boundary. When A = — 1, a new zero mode (and 
eigenvalue 1 of Z) has emerged. 

Note that the eigenvectors in this example have 
no connection to space-time physics and there is 
no association with underlying physical structures, 
like in QCD. However the observed behavior of the 
eigenvalues is quite characteristic. An exact GW- 
operator reacts like a step function. The real eigen- 
value pair disappears and a complex conjugate pair 
of eigenvalues appears. For an approximative GW- 
operator the eigenvalues flow as in the example. It 
is of high interest to quantify the flow and to iden- 
tify, what resolution the respective Dirac operators 
have with regard to the instanton size. 



3. "Artificial" and realistic instantons 

We construct gauge configurations with artifi- 
cial instantons as discussed in jlj, |l^ and study 
for given lattice size 16^ such configurations for in- 
stantons of varying size R (Fig. 2). Real eigenval- 
ues may be associated with instanton modes. 

We recover the typical features of our test ex- 
ample: The Wilson operator reacts quite slowly on 
the change of size and "looses" the instanton at 



Fig. 3. The distribution of the real eigenvalues x of the 
chirally improved operator on 16^ lattices (200 configura- 
tions for each value of the gauge coupling (3). 

size 0(2), the overlap operator looses the instan- 
ton around R ~ 1.5. The chirally improved oper- 
ator is for i? > 3 similar to the overlap operator, 
but identifies also smaller instantons. 

The chirally improved operator may be used ei- 
ther for its own sake, as a Dirac operator, or in or- 
der to improve performance by using it as a start- 
ing point for the determination of the overlap Dirac 
operator. From Fig. 2 one finds, that the improved 
operator is clearly superior to the Wilson operator, 
since it spends less time (fewer configurations) in 
the central region due to its steeper slope. 

In Fig. 3 we show the distribution of eigenvalues 
of the chirally improved operator. The gauge con- 
figurations have been obtained with the Liischer- 
Weisz action [|l6| as detailed in Compared to 
the Wilson action the eigenvalues are much more 
localized around and therefore H will have much 
fewer zero modes. This ought to improve the con- 
vergence properties significantly. However, in our 
computations we found an improvement factor of 
only 2 in the necessary order of the Chebychev 
polynomials. This is likely due to the explicit treat- 
ment of the sector of small eigenvalues of H . Due 
to the numerical complexity of the improved oper- 
ator (it has more coupling terms and is typically a 
factor of 20 more expensive than Wilson's action) , 
we find that using it as input for the overlap opera- 
tor is not economically warranted. However, to use 
it as an approximate Ginsparg- Wilson-type oper- 
ator may be quite successful, as we discuss below. 

The (normalized) gauge-invariant densities 
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Fig. 4. Comparing the inverse participation ratio for the 
overlap operator with the chirally improved operator. For 
both values of the gauge coupling we plot the results of 
200 gauge configurations with a single zero mode. 



(2) 



(implied summation over color and Dirac indices) 
are determined from the zero mode eigenvector 
IV'o); 7a is one of the 16 7-matrices of the 4D- 
Chfford algebra with 70 = 1 . Comparing po and ps 
as determined for a realistic gauge configuration 
with one zero-mode for the three Dirac operators 
considered we find clear instanton bumps at the 
same position in all cases. However, Fig. 2 pointed 
at a varying sensitivity of different Dirac operators 
with regard to identifying instantons. The space- 
time integral /q = / d'^x Pq{x) defines the so-called 
inverse participation ratio and constitutes a mea- 
sure of localization. For instantons with radius R 
one expects Iq oc 

In Fig. 4 we compare Iq for realistic gauge con- 
figurations. We find that the overlap operator 
typically underestimates the localization of the 
zero mode, confirming earlier findings for artificial 
instantons ||l0[] . A possible interpretation lies in 
the "size" of the different lattice Dirac operators. 
The matrix elements of the overlap operator fall 
off like roughly ^ l/2l^^^l whereas the chi- 
rally improved operator is ultra-local (i.e. with 
non- vanishing coupling only for — < \/5)- 
At typical lattice spacings 0(0.1 fm) of today's 
simulations the chirally improved operator ap- 
pears to keep better track of structures smaller 
than 0(0.3 fm). Closer to the continuum limit, at 
smaller lattice spacings, this feature of the overlap 
operator may be no problem, though. 
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